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ABSTRACT 

The massive black hole in our galactic center, Sgr A*, accretes only a small fraction of the gas 
available at its Bondi radius. The physical processes determining this accretion rate remain 
unknown, partly due to a lack of observational constraints on the gas at distances between 
~ 10 and ~ 10^ Schwarzschild radii {Rf) from the black hole. Recent infrared observations 
identify low-mass gas clouds, G1 and G2, moving on highly eccentric, nearly co-planar 
orbits through the accretion flow around Sgr A*. Although it is not yet clear whether these 
objects contain embedded stars, their extended gaseous envelopes evolve independently as 
gas clouds. In this paper we attempt to use these gas clouds to constrain the properties of 
the accretion flow at ~ 10^ R^. Assuming that G1 and G2 follow the same trajectory, we 
model the small differences in their orbital parameters as evolution resulting from interaction 
with the background flow. We find evolution consistent with the G-clouds originating in the 
clockwise disk. Our analysis enables the first unique determination of the rotation axis of the 
accretion flow: we localize the rotation axis to within 20°, finding an orientation consistent 
with the parsec-scale jet identified in x-ray observations and with the circumnuclear disk, 
a massive torus of molecular gas ~ 1.5 pc from Sgr A*. This suggests that the gas in the 
accretion flow comes predominantly from the circumnuclear disk, rather than the winds of 
stars in the young clockwise disk. This result will be tested by the Event Horizon Telescope 
within the next year. Our model also makes testable predictions for the orbital evolution of 
G1 and G2, falsifiable on a 5-10 year timescale. 
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1 INTRODUCTION 

The accretion flows feeding massive black holes have important 
implications for galaxy evolution, for the growth of massive black 
holes, for jet dynamics, and for black hole physics in general. In our 
galactic center, x-ray measurements determine the density and tem¬ 
perature of the gas at the outer edge of the accretion flow (Baganoff 
et al. 2003; Quataert 2002, 2004); the simplest model for time- 
steady, non-rotating, adiabatic accretion (Bondi 1952) then predicts 
an inflow rate Mbondi ~ 10”^ Mg/year. However, this predicted ac¬ 
cretion rate is inconsistent with almost every other measurement 
of the gas in the galactic center, including its bolometric luminos¬ 
ity (Quataert 2002), its rotation measure near the black hole (Agol 
2000; Quataert & Gruzinov 2000; Marrone et al. 2007), and the den¬ 
sity profile inferred from x-ray spectroscopy near the Bondi radius 
(Wang et al. 2013). Extensions to the Bondi model have therefore 
been proposed, including convection (e. g. Quataert & Gruzinov 
2000), rotation (e. g. Narayan & Yi 1995; Narayan & Fabian 2011), 
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magnetic tension (Pen et al. 2003; Pang et al. 2011), and outflows 
(e. g. Blandford & Begelman 1999; Quataert & Gruzinov 2000; 
Yuan et al. 2003). These theories lead to different predictions for 
the density and rotation profiles of the gas in the galactic center; 
however, a lack of observational probes of the gas at intermediate 
radii (between ~ 10 and ~ 10^ Schwarzschild radii, R^) has made it 
difficult to differentiate these models observationally. 

Gillessen et al. (2012) reported the discovery of G2, a gas 
cloud on an extremely eccentric orbit (e ~ 0.98) about Sgr A*, the 
massive black hole in the center of our galaxy. To date, much of 
the theoretical and observational work on G2 has focused on the 
possibility that it will disrupt near pericenter (~ 10^ Rs) and drive an 
episode of increased accretion onto the black hole. Equally exciting, 
however, is the opportunity to use measurements of G2 to study the 
accretion flow feeding Sgr A*: any interaction observed between 
G2 and the ambient gas it moves through amounts to an important 
detection of the galactic center accretion flow. This interaction 
could take the form of a bow-shock ahead of the cloud (Narayan 
et al. 2012; S^dowski et al. 2013a,b; Yusef-Zadeh & Wardle 2013; 
Crumley & Kumar 2013), hydrodynamic disruption of the cloud 
by shear instabilities (Burkert et al. 2012), or deviation from a 
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Keplerian orbit caused by a (magneto-)hydrodynamic drag force 
(Pfuhl et al. 2015; MeCourt et al. 2014). 

Pfuhl et al. (2015) present new observations of Gl, a gas 
cloud similar to G2 which was discovered a decade ago (Clenet 
et al. 2005; Ghez et al. 2005). They find an orbital solution for 
Gl strikingly similar to that of G2; as we discuss in the following 
section, this suggests a common origin for the two objects, with Gl 
preceding G2 by about 13 years. Gl has a slightly lower semi-major 
axis and eecentricity than G2, which Pfuhl et al. (2015) explain as 
the result of a simple drag force due to interaction with the ambient 
gas. If true, this represents a detection of the gas in the galactic 
center, which has so far only been measured at radii much larger, 
and at radii much smaller, than G2’s position. Using this drag force 
to infer the properties of the accretion flow at radii of hundreds 
of AU would provide a valuable test of the various theories for 
the galaetie eenter aeeretion flow. We explore this possibility here; 
in particular, we show that the different components of the drag 
force separately constrain the density and rotation profiles of the 
aeeretion flow. 

Our results will be especially interesting in coming years as 
millimeter wave interferometers resolve event horizon-seale struc¬ 
ture around the black hole (Doeleman et al. 2008; Fish et al. 2011). 
Numerical simulations have shown that the gas dynamics very near 
the black hole is in fact sensitive to the magnetic flux accreted 
from larger radii (e. g. Tchekhovskoy et al. 2011; McKinney et al. 
2012). The gas near G2’s position thus serves as an outer boundary 
condition for the horizon-scale accretion disk, and in many ways 
controls the dynamics that the Event-Horizon Telescope (EHT) 
will observe. Moreover, we show in seetion 3.2 that the drag foree 
inferred from G2’s orbit implies a particular rotation axis for the 
accretion flow. Since this rotation axis is not expected to evolve 
strongly between radii of hundreds of AU and the event horizon,^ 
this is a prediction that will be tested in the next few years when 
the EHT images the accretion disk surrounding Sgr A* (Dexter 
et al. 2010; Broderick et al. 2011; Dexter et al. 2012; Psaltis et al. 
2015; Chan et al. 2015). We discuss this in more detail in section 4. 
Another test will come in early 2018 with the peri-center passage 
of the S2/S0-2 star which may drive a shoek through the aeeretion 
flow visible in x-ray emission (Giannios & Sironi 2013). 

A primary uncertainty in using G2 to probe the galactic center 
accretion flow is that we don’t yet know the nature, origin, or 
structure of the object. G2 is detected as an extended (~ 100 AU) 
cloud in recombination-line emission from hydrogen and helium 
at ~ 10"^ K and as a ~ 1 AU point source in continuum emission, 
interpreted to be thermally emitting dust at ~ 500 K. It is an unusual 
object and its nature and origin are topies of ongoing investigation 
(e. g. Burkert et al. 2012; Murray-Clay & Loeb 2012; Phifer et al. 
2013; Guillochon et al. 2014; De Colle et al. 2014); in particular, 
it’s not yet clear whether the thermal emission represents a star 
embedded within the cloud (cf. Phifer et al. 2013; Witzel et al. 
2014; Valencia-S. et al. 2014). Regardless of its origin, however, 
the ionized component of G2 is spatially extended (~ 100 AU) and 
is measurably distorted by the tidal field of the black hole (Gillessen 
et al. 2012, 2013a,b; Pfuhl et al. 2015). This ionized gas is much 
larger in extent than the tidal radius (~ 1 AU) of any stellar-mass 
companion, so it behaves independently as an extended cloud. In 
this paper, we foeus on the dynamies of the ionized gas alone and 


^ unless the black hole spin is misaligned with the accretion disk (McKinney 
et al. 2013). 


Table 1. (top): Angular momenta for Gl and G2 
inferred from astrometry and velocity data us¬ 
ing the MCMC analysis described in appendix A. 
Units of j are pc^/yr, and error bars indicate 1-cr 
confidence intervals, (bottom): components of the 
torque vector derived from the MCMC chain via 
T ~ (/Gl - iG2)/(tperi, G2 “ tperi,Gl). We also list 
the perpendicular and parallel components of the 
torque, calculated using equations la and lb, re¬ 
spectively. Torques are normalized to the angular 
momenta so the units are [r/ j] = yr“k As above, 
error bars indicate the 1-cr confidence interval. 
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ignore the central thermal emission. We thus remain agnostic about 
the formation scenario for G2. 

We have organized this paper as follows: in section 2 we show 
that the measured orientations of Gl and G2 imply a drag force with 
a significant component out of the orbital plane; such a drag force 
can only be caused if the background G2 moves through rotates at a 
significant fraction of the local Keplerian speed. We present a more 
detailed analysis in section 3, where we write down a simplified 
equation of motion for G2 and numerically integrate its trajectory. 
We show that we can reproduce the observed measurements for 
Gl and G2 with reasonable choices of our model parameters. We 
conclude in section 4 by discussing the possibility of using future, 
improved measurements of Gl and G2 to further constrain the 
properties of the accretion flow in the galactic center. We also 
discuss how our model can be tested using upcoming measurements 
by the Event Horizon Telescope. 


2 THE TORQUE CONNECTING Gl AND G2 

Our analysis is predicated on the assumption that Gl and G2 follow 
the same trajectory, as suggested by Pfuhl et al. (2015). We begin 
by discussing why we expect this is a reasonable assumption. Gl 
and G2 have strikingly similar properties and are on unusual, highly 
eccentric orbits (Pfuhl et al. 2015). The angular momenta of Gl and 
G2 are aligned to within (16 + 3)°; if the orbits were isotropically 
distributed on the sky, there is only a ~ 2% probability for such 
an alignment to occur by chance. Compounding this probability is 
the fact that the eccentricity vectors of Gl and G2 are also aligned 
to within (29 + 4)°; accounting for the fact that the eccentricity 
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Figure 1. Visualization of the torque vectors presented in table 1. His¬ 
tograms show the probability distributions for the x- (red), y- (green), and 
z- (blue) components of the torque vectors computed from the MCMC fit 
to the data in appendix A. The data constrain the perpendicular component 
of the torque t±_ / j somewhat more strongly then r alone, but there is a 
reasonably strong detection of a non-zero torque in both cases. This torque 
is oriented in a way that both re-orients and increases the angular momen¬ 
tum of the cloud; such a torque cannot be produced by a drag force acting 
on a spherical object moving through a stationary background. However, 
we show in section 3 that it can be caused by motion through a rotating 
accretion flow. 


and angular momentum vectors are orthogonal, there is a ~ 0.04% 
chance of drawing two such aligned orbits from an isotropic dis¬ 
tribution. (This probability remains small, ~ 1%, even if G1 and 
G2 instead originate from random locations in the clockwise stellar 
disk.) Moreover, out of the ~ 400 year orbital period of G2, the 
pericenter passages of G1 and G2 coincide to within ~ 13 years; 
such temporal alignment has only a < 5% probability of occurring 
at random. Together, these probabilities imply that, if G1 and G2 
are unrelated objects on independent orbits, we should expect to see 
> 10"^ similar objects within ~ 0.04pc of Sgr A*. Such a population 
would be readily detectable, even accounting for possible selection 
biases for clouds with large line-of-sight velocities relative to earth 
and with epochs of pericenter passage near the present day. 

An alternative scenario in which G1 and G2 are related, but 
somehow following different orbits, also seems improbable. Since 
G1 and G2 are on such eccentric orbits, they spend most of their 
time close to apocenter and likely formed there. But, if G1 and 
G2 follow Keplerian orbits, their apocenter locations differ signif¬ 
icantly due to their high eccentricities and differing orientations 
and semi-major axes (see, e. g. figure 4). Their orbital periods also 
differ by a factor of ~2-3. So if G1 and G2 follow separate orbits, 
they must have formed in very different places and at very differ¬ 
ent times, only to come into co-location and alignment near their 
pericenter passages; such a scenario seems contrived. (We note that 



z (pc) 


Figure 2. Diagram sketching G2’s orbit, along with the perpendicular torque 
vector required to rotate its orbital plane into Gl’s orientation. Gray arrows 
indicate the 1-cr uncertainty in the direction of the torque vector. This torque 
implies a force out of the page in the upper-left quadrant of the figure (post- 
pericenter) and into the page in the lower quadrants (pre-pericenter). A 
rotating accretion flow with its rotation axis tilted in the direction of can 
naturally explain this out-of-plane drag force, along with the reversal of the 
force near pericenter. 


the scenario proposed by Guillochon et al. 2014, in which G1 and 
G2 formed near pericenter, does not suffer from this problem; our 
analysis is consistent with that formation scenario.) 

Hence, in this paper we assume that G1 and G2 follow the 
same trajectory, with G1 preceding G2 by ~ 13 years (see Pfuhl 
et al. 2015). If traveling through a vacuum, this trajectory would be 
Keplerian. However, the trajectory evolves due to interaction with 
ambient gas in the accretion flow. The small differences in their 
orbital parameters thus encode information about the drag force 
exerted on the trajectory and properties of the accretion flow. We 
model this evolution and attempt to use it to constrain the properties 
of the accretion flow. 

The orientations of the orbital planes for G2 and G1 are known 
somewhat more precisely than the energies of their orbits (Gillessen 
et al. 2012, 2013b; Phifer et al. 2013; Pfuhl et al. 2015); this con¬ 
strains the direction of the angular momentum vector better than its 
magnitude. Hence, we separately consider two components of the 
torque vector r: 


/ = [g-j «>}•]• I= y (la) 

/ = j ■ / = T||, (lb) 

where g is the metric tensor (the identity matrix in Cartesian coordi¬ 
nates), j is a unit vector in the direction of the angular momentum 
7 , and j is its magnitude. The perpendicular torque reorients 
the angular momentum j\ while the parallel torque Ty changes it 
magnitude. 

We estimate the angular momentum vectors for G2 and 
G1 from the astrometry and velocity measurements published in 
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Gillessen et al. (2013a,b)^ and in Pfuhl et al. (2015). We then cal- 
culate the components of the torque vector ry and / j implied 
by the change in angular momentum. We provide details of this 
fit in appendix A. Table 1 lists the angular momentum and torque 
vectors we obtain. Figure 1 visualizes this information to illustrate 
the uncertainty in different components of r and . Though the 
uncertainties are large, there is good evidence for a change in orbital 
angular momentum between G2 and Gl. 

Both the perpendicular component, r^/ j, and the parallel 
component, t\\/ 7 , of the torque provide evidence for a rotating ac¬ 
cretion flow. We first discuss the perpendicular component, which 
re-orients the orbit. Physically, this torque corresponds to the com¬ 
ponent of the drag force out of G2’s orbital plane. A symmetric 
object moving through a stationary background cannot experience 
such a force; our finding that |t^| ~ r thus implies either that G2 
is flattened like an airplane wing to generate lift, or that it moves 
through a rotating accretion flow with a rotation velocity compa¬ 
rable to G2’s orbital velocity; in fact, we show in section 3.2 that 
both conditions are necessary to fit the data. If G2’s post-pericenter 
evolution indeed follows Gl, this is a detection of rotation in the 
galactic center accretion flow. Figure 2 illustrates the geometry of 
the inferred drag force. 

The parallel component of the torque Ty/ j, which corresponds 
to the component of the drag force in G2’s orbital plane, also 
implies a rotating accretion flow in the galactic center. For the 
orbital parameters quoted in Gillessen et al. (2013b) and in Pfuhl 
et al. (2015), Gl has a larger angular momentum than G2. If a drag 
force can turn one orbit into the other, it must spin up the orbital 
angular momentum of the cloud (ry > 0 ), rather than torquing 
it down; this is possible only if the background rotates rapidly 
with a velocity comparable to G2’s orbital velocity. Pfuhl et al. 
(2015) neglected rotation, and thus could not fit the increase in 
angular momentum between G2 and Gl: they found a larger angular 
momentum for G2, and smaller energy and angular momentum for 
Gl, than are observed. If the published orbital elements for Gl and 
G2 prove accurate, and if G2 continues to follow Gl’s orbit over the 
next ~ 10 years, this is again a detection of rotation in the galactic 
center. 

The inferred torque vector shown in table 1 and in figure 1 con¬ 
stitutes the essential result underpinning our analysis. We quantify 
it in the next section using numerical integrations and an idealized 
model for G2’s evolution. 


3 NUMERICAL INTEGRATIONS 
3.1 Method 

We model the gas cloud as a point particle moving in the Keplerian 
potential of the black hole, subject only to a drag force due to its 
interaction with ambient gas in the accretion flow. We ignore the 
stellar contribution to the potential since the timescales for orbital 
precession and gravitational relaxation are much longer than the 
orbital period of G2 (Kocsis & Tremaine 2011). 

We adopt a drag force oc + 2/[j3M^]) as pre¬ 
sented in McCourt et al. (2014), where p = quantifies 

the magnetic field strength in the background plasma, Pbg(^) is the 
density of the background gas, Vrei = v - Vbg(r) is the relative veloc¬ 
ity between the cloud and the background gas, and M is the Mach 


^ made publicly available at https://wiki.mpe.mpg.de/gascloud/ 
PlotsNData 


number. Physically, the last term represents an enhancement over 
the hydrodynamic drag due to the magnetic field in the background 
plasma. This leads to the following approximate equation of motion 
for the cloud: 


d^r ^ GM.r Pbg(r) / ^ 2 \ 

df r3 Meioudl PM^j 

The ram pressure vector R in equation 2 is 

^ ~ '^cloudi'^rel ’ V) V + ^cloud-^cloud(Prel ’ ^) d 
"l"'^cloud-^cloud(Prel ’ ^) d. 


( 2 ) 


( 3 ) 


where v, d, and b form a right-handed coordinate system aligned 
with the cloud. We thus allow the cloud to be elongated along its 
orbit with a (constant) aspect ratio Ldoud/^cioud- The ram pressure 
vector R reduces to ^cioud^rei^rci when v^i is parallel to v. If the 
cloud is quasi-spherical, with Ldoud = ^cioud. the drag force points 
approximately in the direction opposite the relative velocity Vj-d*^ A 
non-spherical cloud develops a drag force in a direction different 
from Vrd; this force is akin to aerodynamic lift and we show in 
section 3.2 that it is essential to fit the data. 

Throughout our analysis, we adopt the best-fit values for the 
mass of the black hole, M, = 4.3 x 10^ M©, and for the distance 
to the galactic center, Rq = 8.3 kpc, from Gillessen et al. (2009). 
We also take Mdoud ~ 3 Mearth and Rdoud ~ 100 AU, as suggested 
in Gillessen et al. (2012). We do not vary these parameters when 
attempting to fit our model to the data. 

In order to derive an initial position and velocity, we took the 
astrometry and velocity data for G2 from Gillessen et al. (2013a,b). 
We fit the observed positions and line-of-sight velocities with 
second-order polynomials in time and use these fits to produce 
an initial condition for G2 at 2013.33. In units of parsecs and years, 
this reads: 


r (2013.33) = 

V (2013.33) = 


0.003537 

-0.00002680 

-0.001088 

-0.002075 

0.0008509 

0.002189 


(4a) 


(4b) 


where the origin is centered on Sgr A* and the v-, y-, and z-axes 
correspond to right ascension, declination, and the line-of-sight 
distance, respectively. This initial condition is consistent with the 
published Kepler elements for G2’s orbit (Gillessen et al. 2012; 
Phifer et al. 2013). Starting from this initial condition, wc integrate 
the equations of motion using the NDSolve function in Mathemat- 
ica. We have made our code publicly available on github."^ 

The background accretion flow enters into our equation of mo¬ 
tion through the density Phg(r), the relative velocity v^i = v - Vbg(r), 
the magnetic field strength p, and the Mach number M, which de¬ 
pends on the temperature of the background gas. Unfortunately, 
there are very few observational constraints on the background 
accretion flow at the range of radii where G2 orbits. We proceed 
by adopting a model for the background gas; we then compare 
resulting trajectories with measurements of Gl and G2 to constrain 
the parameters in our model. In this paper, we adopt an extremely 


^ not precisely opposite, because equation 3 is appropriate for a cloud with 
flat surfaces, rather than a sphere. Due to the lack of perfect symmetry, 
non-spherical clouds can experience lift. 

^ https://github.com/mkmcc/g2-drag-force 
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simplistic model for the density, temperature, and velocity of the 
baekground gas: 


lrV“ 

Pbg(r) = Po - 
\ro 

(5a) 

GM. 

rbg(r) = — 

(5b) 


VbgW-Zkep^ r j 7r ’ 

(5c) 


where po = 1.7 x 10^ ^Earth/pc^ and fq = 0.04 pe are fixed by 
x-ray measurements of Bremsstrahlung emission near the Bondi 
radius (Quataert 2002, 2004). This model is specified by the power- 
law exponent of, by the rotation parameter /kep = Vbg/Vcirc, by the 
two angular coordinates (6, 0) of the angular momentum vector 
J oc [sin(^)cos(0), sin(^) sin(0), cos(^)], and additionally by the 
assumed magnetie field strength (via p in equation 2). This model 
does not represent an internally eonsistent theory for the aeeretion 
flow in the galactic center (see, e. g.. Yuan et al. 2003); rather, it is a 
minimal model which might rotate the orbital plane of an infalling 
gas eloud. It is straightforward to extend our analysis by replaeing 
equation 5 with a suitable model for the gas in the galaetie eenter; 
the analysis discussed in the next section might then be used to 
differentiate among the various theories for the accretion flow. Due 
to the limited amount of data for G1 and G2 and the uncertainty 
in their masses and shapes, we cannot yet do so. However, we 
describe in section 4 how continued observations of G1 and G2 
eould improve our eonstraints on a timeseale of ~5-10 years. 

After integrating each trajectory, we estimate its likelihood 
assuming the probability P oc Qxp(-x^), where 

9 [model(t/) - data/]^ 

is ealeulated using the published astrometry and veloeity measure¬ 
ments for G1 (Pfuhl et al. 2015) and G2 (Gillessen et al. 2013a,b). 
We maximize this likelihood over all initial conditions and back¬ 
ground models. We use two complimentary methods for this opti¬ 
mization. First, we use Mathematica’s built-in simulated annealing 
optimizer, whieh produces a population of orbits tightly clustered 
around the best-fitting solution; we show this best-fitting orbit in 
figures 3, 5, and 7. Second, we use the emcee MCMC optimizer 
(Foreman-Maekey et al. 2013), whieh is an implementation of the 
affine-invariant algorithm presented in Goodman & Weare (2010). 
We use the resulting Markov chain to estimate the uncertainty in 
our model parameters (figures 6 and 7); however we caution that 
the uncertainty in our results is likely dominated by simplifica¬ 
tions in our model, not by the statistical noise quantified by the 
Markov chain. By analogy with the Gaussian distribution, we define 
the “1-cr” uncertainty for the parameters in our Markov chain as 
the interval between the d and \ - d quantiles of the data, where 
d = ierfc(l/V2) ~ 0.159. For the “2 -cr” uncertainty interval, 

d - ^erfe( V 2 ) ~ 0.0228. The best-fit values we quote for the 
Markov chain are the most probable values determined via the 
binning method deseribed by Knuth (2006). 

We find the best-fitting orbit has ;^f^/d.o.f. ~ 4, rather than 
a value ~ 1 expected for a good fit to the data. As discussed in 
appendix A and in the following seetion, this is eomparable to the 
1^.0. f. we obtain by fitting G1 and G2 to separate Keplerian or¬ 
bits; it thus seems likely that the errors quoted for the observational 
measurements are underestimated by a factor of < 2. Following 
Guillochon et al. (2014), we quantify this additional error by re¬ 


peating our analysis with the modified statistie: 

[model(t/) - data/]^ 


Pocfl-^expl- 

i J27Tsf 


(7) 


where 


= cr^ + /^(data)^ 


( 8 ) 


and the extra parameter / represents a possible systematic uncer¬ 
tainty in addition to the error-bar cr quoted for the observational 
points. We allow different /’s for astrometry and line-of-sight veloc¬ 
ity measurements. This statistic produces similar results to the;^^^ 
statistie mentioned above, but yields a more reasonable likelihood 
with an equivalent;^^/d.o.f. ~ 1.5. 

In our optimization, we assume flat priors on a, /kep, and 
on perturbations to the initial condition. We assume the Jeffreys 
prior for the parameter /. The polar and azimuthal angles 6 and 
0 define the orientation of the rotation axis; we assume these are 
isotropically distributed, i. e. with flat priors for sin(^) and for 0. We 
also assume flat priors for log p and for log (Ldoud/^cioud)- While 
these are not neeessarily “uninformative” priors, different choiees 
do not appear to bias our results by more than ~ 1-cr. Different 
choices for the background model in equation 5 have a larger effect. 

Equations 2-5 fully specify our numerical model for G2’s 
evolution, and equations 6-8 specify our comparison to the data. 
We show our results in the following seetion. 


3.2 Results 

The dark blue curve in figure 3 shows the best-fitting orbit we found 
with the drag force model in equation 2. In this figure, the red points 
indieate the Br-y observations of G2 from Gillessen et al. (2013b) 
and the green and blue points respectively show the L-band and 
Br-y observations of G1 from Pfuhl et al. (2015). This solution 
has a;^^/d.o.f. = {4.2,1.4}, where the first number is a standard 
defined using equation 6 and the second number is calculated using 
the larger error bars found via equation 7. (We show these larger 
error bars in light colors in figure 3.) This quality of fit is essentially 
equivalent to the best fit we obtained while fitting G1 and G2 to 
separate Kepler orbits (y^/d.o.f. = (4.4,1.5); see appendix A); the 
data thus do not presently favor our model over one in which G1 
and G2 are unrelated objects on distinct Keplerian orbits. However, 
as discussed in section 2, G1 and G2 are on unusual orbits and have 
very different apocenter locations; if they are unrelated objects, they 
formed in different plaees and at different times, and it thus seems 
highly improbable to find them so elosely eo-loeated at perieenter. 
This low probability is not aeeounted for in the;^^ statistie quoted 
above, but should be considered when evaluating a model in which 
G1 and G2 are unrelated. 

The gray curve in figure 3 shows a model in which the drag 
force is restricted to the orbital plane, as in Pfuhl et al. (2015). 
This model fits the data similarly to ours (y^/d.o.f. = (5.6,1.7|), 
but allowing for rotation more elosely matehes the post-perieenter 
radial velocity. While these two models both fit the existing data, 
they predict very different evolution following the first perieenter 
passage. For Gl, which is ahead of G2 by about a deeade, the two 
models will begin to diverge appreeiably in 5-10 years. 

Figure 4 shows a perspective rendering of our best-fit orbit, 
along with the Kepler ellipses fit to Gl and G2 by Pfuhl et al. (2015) 
and Gillessen et al. (2013b). While the cloud’s trajectory does lie 
tangent to the Kepler ellipses for Gl and G2, this alignment is 
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Figure 3. (left:) Comparison of best-fit model (blue line) with astrometric data. Smaller, dark error bars correspond to those reported in Gillessen et al. 
(2013a,b) and Pfuhl et al. (2015). Larger error bars include the systematic error found from our maximum-likelihood analysis (see section 3.1). Though the 
cloud appears to be swept away in this figure, this is a projection effect; the orbit does shrink and circularize (see figure 5). The light gray orbit shows a model 
in which the drag force is restricted to the plane of the orbit, as presented in Pfuhl et al. (2015). The location of Sgr A* is marked with an asterisk, (right:) 
Comparison of best-fit model with line-of-sight velocity. As in the astrometry plot, the smaller error bars are from Gillessen et al. (2013a,b) and Pfuhl et al. 
(2015), while the larger ones account for the systematic error determined by our maximum-likelihood analysis. In both panels, the filled black circles mark the 
location of G1 in 2018 around the time of the peri-center passage of S2/S0-2, for the Pfuhl et al. (2015) model. The differences with our model (blue line) will 
be easily detectable. 



Figure 4. Perspective rendering of the best-fit orbit. Gray ellipses show the 
Keplerian fits for G2 and for Gl. In our model, the gas cloud initially comes 
in from apocenter approximately following the Keplerian orbit fit to G2. 
The cloud’s orbital plane is misaligned with the rotation of the background 
flow, resulting in a perpendicular component to the drag force. The resulting 
torque rapidly re-orients the orbital plane and, after ~ 5 pericenter passages, 
the orbital plane is aligned with the rotation axis of the accretion flow. No 
further reorientation occurs, but the orbit circularizes and the cloud becomes 
approximately co-moving with the rotating gas. The drag force thus drops 
significantly; the cloud slowly sinks in toward the black hole. Figure 5 
shows this orbit from several different angles to better illustrate the rotation. 


transitory. This figure is drawn in a coordinate system with decli¬ 
nation increasing “up” the page, line-of-sight distance increasing 
to the right, and right ascension increasing into the page. Figure 5 
shows the same orbit in a rotated coordinate system aligned such 


that the cloud is initially on the x-axis, and the rotation axis of 
the background flow lies in the x-z plane. This coordinate system 
more clearly shows the rotation of the cloud’s orbital plane. In 
these figures, the gas cloud initially comes in from apocenter on 
a nearly radial orbit. The cloud’s orbital plane is misaligned with 
the rotation of the background flow, resulting in a perpendicular 
component to the drag force. As the cloud approaches the black 
hole, this perpendicular drag force increases due to both the increas¬ 
ing background density and the increasing rotation velocity. The 
resulting torque rapidly re-orients the orbital plane and, after ~ 5 
pericenter passages, the orbital plane is aligned with the rotation 
axis of the accretion flow. No further reorientation occurs, but the 
orbit circularizes and the cloud becomes approximately co-moving 
with the rotating gas. The drag force thus drops significantly, and 
the cloud slowly spirals in toward the black hole. Of course, after 
many pericenter passages, the cloud may be tidally distorted to 
the point where it intersects itself; this process could lead to much 
more rapid circularization and inflow, but is not accounted for in 
our calculation. 

Figure 6 shows probability distributions of the parameters in 
our model, assuming the probability defined in equation 7 (the 
standard analysis yielded similar results, but with somewhat 
smaller uncertainties). The data constrain some parameters much 
more strongly than others: for example, we localize the orientation 
of the rotation axis of the accretion flow to within ~ 20°. The rota¬ 
tion axis for the accretion flow we find corresponds to a polar angle 
^ ~ (111 + 16)° and an azimuthal angle 0 ~ (-32 + 23)°. This is 
misaligned with the stellar disk by almost 90°, and with G2’s orbit 
by about 60°. Our analysis does not yet provide an easy identifica- 
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Figure 5. Perspective rendering of the best-fit orbit. In order to more clearly show the reorientation of the orbit, the coordinate system is rotated so that 
apocenter lies along the x-axis and the rotation axis of the accretion flow lies in the x-z plane. The cloud initially plunges toward the black hole, and the drag 
force quickly re-orients the orbital plane after ~5 pericenter passages. Once the orbit aligns with the accretion flow, it becomes approximately co-moving with 
the gas. The drag force is thus much reduced, and the orbit circularizes and slowly moves in from that radius. 



Figure 6. Probability distributions of the parameters in our model, as determined by our maximum-likelihood analysis (we obtained very similar results from 
the;^'^ analysis described in section 3.1). The data constrain some parameters much more strongly than others: for example, we localize the orientation of the 
rotation axis of the accretion flow to within ~ 20° {left panel). The parameters a and both control the magnitude of the drag force; we thus And they are 
strongly degenerate {center panel). The magnetic field strength can be determined by other observations, however (e. g. Eatough et al. 2013); this would enable 
a constraint on the gas density in the galactic center. Similarly, the rotation parameter /^ep and the cloud aspect ratio Tdoud/^cioud control the component of the 
drag force out of the plane; these parameters are also strongly degenerate {right panel). The shape of the cloud can be inferred from modelling its velocity 
shear as determined by Br-y observations; this would enable a constraint on the rotation profile of the accretion flow. 


tion for the origin of rotation in the galactic center. Identifying the 
source of rotation could be further complicated by the back-reaction 
of the gas cloud onto the accretion flow: since the combined mass 
of G1 and G2 is comparable to the total mass expected to be in the 
accretion flow, the drag force on G2’s orbit could have rotated the 
accretion flow by several degrees from its original position. How 
the flow responds to this back-reaction depends in detail on mo¬ 
mentum transport within the gas; viscosity, convection, turbulence, 
and magnetic fields could all influence the process significantly. We 
have not modeled any such back-reaction in this paper; however it 
would be straightforward to add once a theory for the background 
flow is specified. 

While our model constrains the rotation axis of the accretion 
flow, figure 6 shows that other combinations of parameters are 
degenerate. The density of the background gas (governed by its 
power-law slope a) and the magnetic field strength (quantified by 
J3) both affect the magnitude of the drag force. We can increase the 
drag force by either raising the density of the gas, so that the cloud 
accelerates a denser column of ambient gas, or by increasing the 
field strength, so that magnetic tension couples the cloud to a larger 
volume of background gas (McCourt et al. 2014). Since the evolu¬ 


tion of the orbit depends only on the overall magnitude of the drag 
force, we find that the parameters a and p are almost completely 
degenerate. Similarly, the rotation parameter /kep and the cloud 
aspect ratio Ldoud/^cioud independently control the component of 
the drag force out of the orbital plane; this perpendicular drag force 
is proportional to x Ldoud/^cioud- Therefore, we only constrain 
their product and these parameters are also strongly degenerate. 

Thus our model does not directly constrain cr, p, /kep, or the 
shape of the cloud. However, some of these parameters may be 
constrained by complementary observations. The shape of the cloud 
can be inferred from Br-y observations of the velocity shear; the 
measurements in Gillessen et al. (2013b) suggest Ldoud/^cioud ~ 10- 
20. This would imply a fairly large rotation rate; /kep ~ 0.5-1.0. 
Similarly, observations of the magnetar in the galactic center (Ken- 
nea et al. 2013) suggest a magnetic field strength with ~ 1 at 
large radii (~ 0.1 pc; Eatough et al. 2013). Since p is not expected 
to evolve strongly with radius in many models of the galactic center, 
this suggests a fairly steep density profile with a ~ 0.7-0.9. This 
is within the constraint of the x-ray luminosity, roughly consistent 
with a radiatively inefficient accretion flow (Yuan et al. 2003) and 
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Figure 7 - continued Comparison of the evolution of Kepler elements with 
time for some well-fitting orbits with derived values for G1 and G2. Red 
lines show a random sample of well-fitting solutions. Yellow and blue bands 
indicate 1- and 2-cr intervals. Black points are from Gillessen et al. (2013b) 
and from Pfuhl et al. (2015); we use the data from these papers in our 
analysis. For comparison, the gray points show results from Gillessen et al. 
(2012), Gillessen et al. (2013a), and Phifer et al. (2013). 


with the density profile inferred from x-ray spectroscopy at larger 
radii (Wang et al. 2013). 

Interestingly, if we require G2 to be spherical in shape, we 
cannot reproduce the inferred shift in its orbital plane. Even with 
rotation at the Keplerian speed, matching the out-of-plane drag 
force yields an in-plane force so large that the cloud rapidly spirals 
into the black hole. We can only match the measurements if we 
allow the gas cloud to have a non-spherical shape, so that the drag 
force points in a direction different than the cloud’s relative velocity. 
A misalignment between the drag force and the relative velocity 
is analogous to aerodynamic lift; in a sense, G2 thus partially flies 
through the accretion flow as it falls toward the black hole. The drag 
force we measure implies a glide ratio ~ 2 for G2, comparable to 
that of the space shuttle (~ 4). Our best-fit solutions imply an aspect 
ratio Lcioud/^cioud ~ 10, similar to what is seen in simulations (e. g. 
Burkert et al. 2012; Guillochon et al. 2014) and to what is inferred 
observationally (Gillessen et al. 2013b; Pfuhl et al. 2015). More 
precise determinations of the shapes of G1 and G2 from fitting 
their measured velocity shears could thus constrain the rotation 
parameter /kep; this in turn would provide information about the 
(unknown) thermodynamics of the accretion flow. 

Figure 7 compares the evolution in Kepler elements of well- 
fitting models as a function of time with derived values for G1 and 
G2. The yellow and blue bands indicate the 1- and 2-cr intervals 
derived from the analysis shown in figure 6. (Within the context 
of the model defined in equation 5, the error bars show the range 
allowed by the data. The predictions could of course change by sev¬ 
eral cr under a different evolutionary model.) The red curves show 
a random sample of well-fitting solutions - these yield essentially 
the same and cannot be differentiated using the existing data. 
Figure 7 shows that the energy loss is strongly episodic - the semi¬ 
major axis a decreases sharply at pericenter passages, but is nearly 
constant between them. The angular momentum increases more 
continuously, however: e decreases smoothly with time, as do the 
angles Q! and oj. The inclination /' deceases sharply after the first 
pericenter passage, then remains constant for several decades. This 
provides an alternative test of our model, in addition to the astrom¬ 
etry and velocity curves shown in figure 3. Continued observations 
of e and Q' thus represent a promising way to test this model in the 
coming years and to improve constraints on the parameters for the 
galactic center accretion flow. 

Thus, we have found a model trajectory for gas clouds in the 
galactic center which matches the observed positions and velocities 
of G1 and G2, and which is consistent with available constraints 
on the gas density and magnetic field strength in the galactic center. 
This model implies a particular rotation axis for the accretion flow 
feeding the black hole Sgr A*; in the following section, we describe 
the orientation of this axis relative to other features in the galactic 
center. 
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Figure 8. Comparison of our constraint on the rotation axis (left panel 
of figure 6) with other measurements in the galactic center. The blue and 
orange error bars show constraints on the orientation of the horizon-scale 
accretion disk from measurements by the Event Horizon Telescope (Dexter 
et al. 2010 and Broderick et al. 2011, respectively). Currently, the EHT 
limits the orientation of the rotation axis to one of four possibilities; our 
constraint is consistent with one of them. Red error bars mark the orientation 
and width of the clockwise stellar disk in the galactic center (Bartko et al. 
2009, large error bars, and Yelda et al. 2014, small error bars). Text markers 
indicate the orientations of G2 and G1 (table Al). Though we did not 
include any information about the clockwise stellar disk in our fit, we find 
evolution consistent with the G-clouds originating in that disk. The rotation 
axis we identify lines up closely with the molecular ring in the galactic 
center (green error bars; Jackson et al. 1993); this supports a scenario in 
which the accretion flow originates with gas inflowing from the molecular 
ring. 


4 DISCUSSION 

The dynamics of the accretion flow around the massive black hole 
in our galactic center is an important topic of ongoing research. 
One fundamental problem in this area is that the Bondi accretion 
rate implied by x-ray observations at large radii is much higher 
than the constraints implied by the bolometric luminosity of the 
gas and by its rotation measure (see Yuan & Narayan 2014 for a 
recent review). Most of the gas which could accrete from large 
radii therefore never reaches the black hole. A number of theories 
have been proposed to explain this deflcit, but we can’t currently 
differentiate among them due to a lack of observational constraints 
on the gas in the galactic center. 

Gillessen et al. (2012) reported the discovery of G2, a gas 
cloud on a highly eccentric orbit about the massive black hole in 
our galactic center. G2 has passed through a range of radii where 
differing models for the galactic center accretion flow yield different 
predictions for the density and rotation rate of the gas. Detection 
of any interaction between G2 and this background gas would 
therefore provide a critical test of the various theories of galactic 
center accretion flow. 

Following Pfuhl et al. (2015), we work with the assumption 
that G2 and the similar cloud G1 follow the same trajectory, with 


G1 preceding G2 by ~13 years; we quantify in section 2 why this 
is a reasonable assumption. Though the orbital parameters of G1 
and G2 are strikingly similar, they contain important differences. 
Under the assumption that G2 follows Gl, these small differences 
represent evolution due to interaction with the accretion flow, and in 
turn implies a particular orientation for the drag force acting on the 
cloud (see table 1 or figure 1). This drag force is inconsistent with 
motion through a static medium, and thus implies that the accretion 
flow in the galactic center rotates. We model this drag force in 
more detail in section 3, where we show that we can reproduce 
the observed trajectories of Gl and G2 with a simple model and 
reasonable choices for the properties of the galactic center accretion 
flow (figure 3). 

Figure 6 shows our preliminary results, in which we localize 
the rotation axis of the accretion flow to within 20°. We show this 
result in more detail in figure 8 and over-plot earlier constraints 
on the rotation axis determined by the Event Horizon Telescope 
(EHT; Dexter et al. 2010; Broderick et al. 2011; Psaltis et al. 2015). 
Degeneracies currently limit the EHT determination to one of four 
possibilities; our constraint is consistent with only one of them.^ To 
our knowledge, this paper represents the first unique determination 
of the rotation axis at intermediate radii (~ 10 < ^ < in 

the galactic center. 

In figure 8, we also mark the position of the (inner edge of 
the) clockwise disk with red error bars (Bartko et al. 2009; Lu 
et al. 2009). We note that while we did not include information 
about the clockwise disk in modeling the trajectory of G2, we And 
evolution entirely consistent with an origin within the clockwise 
disk. Yellow arrows sketch the evolution of the clouds in our model: 
they originate in orbits consistent with the clockwise disk and end 
up aligned with the rotation axis of the accretion flow. 

The rotation axis we identify is in close agreement with the 
orientation of G359.944-0.052, a linear x-ray filament which points 
towards Sgr A*. Li et al. (2013) propose this feature results from 
a parsec-scale jet emanating from Sgr A*: the jet collides with 
the Eastern arm of the Sgr A West HII region, driving a shock 
front which is observable in infrared and radio, and accelerating 
relativistic electrons which are seen in x-ray. The long axis of the 
hypothetical jet path corresponds to an azimuthal angle 0 ~ (-35 + 
2)°. If the jet propagates along the rotation axis, this measurement 
is consistent with our finding that 0 ~ (-32 + 23)°. The rotation 
axis is also similar that inferred by jet models of Markoff et al. 
(2007), ^>75°,0~ -25°, and model-dependent constraints on 
the spin axis by Meyer et al. (2007), 0 ~ - 25°. Surprisingly, this 
rotation axis is also nearly coincident with the galaxy’s rotation 
axis, 0 = -31.7° (Reid & Brunthaler 2004). 

Since we find a rotation axis misaligned with both the stellar 
clockwise disk (by ~ 90°), and with G2 (by ~ 60°), our model ar¬ 
gues against scenarios in which gas in the accretion flow originates 
from either of these sources. The rotation axis we identify does 
align quite closely with the circumnuclear disk, a massive torus of 
molecular gas with an inner edge ~ 1.5 pc from Sgr A* (0 = 118° 
+ 15°, 0 = -25° +15°; Jackson et al. 1993). Our analysis suggests 
that the gas in the accretion flow hails from the circumnuclear disk, 
rather than from the young clockwise stellar disk as is commonly 
assumed. 

In addition to constraining the rotation axis of the accretion 
flow, we find degenerate relationships between the gas density and 


^ The angle 0 we use in this paper is related to the angle which is 
measured in degrees east of north, by ^ = 90° - 0. 
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magnetic field strength, and between the rotation profile and the 
shape of the G2 cloud. These degeneracies could be broken with 
further observations, yielding first measurements of the gas density 
and rotation speed at radii of 100s of AU in the galactic center. 
This is an exciting possibility, because these measurements would 
constrain the thermodynamics and momentum transport within the 
accretion flow, providing important information about the physics 
of black hole accretion in low-luminosity systems such as Sgr A*. 

We find that G2’s motion is with the accretion flow rather than 
against it (figure 8; c.f. S^dowski et al. 2013b). With a relative 
inclination of ~60° between G2’s pre-pericenter orbit and the ac¬ 
cretion flow, this is similar to the M3 simulation of Abarca et al. 
(2014). This lowers the relative velocity somewhat, and may help 
explain non-detections of the bow-shock ahead of the cloud (e.g. 
Tsuboi et al. 2015; Bower et al. 2015). A potentially important 
effect we have not included in our analysis is the back-reaction of 
the drag force onto the accretion flow; since the combined mass of 
G1 and G2 is comparable to the total amount of gas in the accre¬ 
tion flow, the drag force acting on the clouds could also influence 
the evolution of the accretion flow. This is a complex process, as 
the back-reaction depends on angular momentum transport in the 
accretion flow and could be affected by magnetic fields, turbulent 
transport, and viscosity; we have therefore chosen not to model it 
in our preliminary study. 

Though our model requires the gas cloud to be tidally elon¬ 
gated in its direction of motion, we have not considered the possible 
disruption of the cloud by tidal forces. While this disruption is a ro¬ 
bust outcome in hydro simulations (e. g. Guillochon et al. 2014), G1 
evidently survived for at least six years after its pericenter passage. 
We therefore assume that some additional process (e. g. magnetic 
tension) enables the clouds to survive pericenter passage. We will in¬ 
vestigate this possibility in the future using magneto-hydrodynamic 
simulations. 

We have not yet considered evolution in the size, shape and 
mass of the cloud. It is straightforward to generalize our model 
and let these values evolve with time, given a theoretical or ob¬ 
servational prescription. However this is not likely to change our 
constraints. The change in orbital parameters of the cloud occurs 
mostly at pericenter where the drag force is strongest. Therefore 
the properties of the cloud far from pericenter have a limited in¬ 
fluence on its orbital evolution. Furthermore, the rotation axis of 
the accretion flow is insensitive to changes in the properties of the 
cloud. 

Our model is not unique; the data for G1 and G2 can be fit 
to independent Keplerian orbits, or by the simpler drag model 
proposed in Pfuhl et al. (2015). In both cases, we argue that our 
model is physically more reasonable. Crucially, these three models 
predict very different evolution for G1 in the next ~5-10 years 
(see figures 3 and 7). Continued observations of G1 are therefore 
essential to test our model and to provide stronger constraints on 
the properties of the accretion flow. If our model proves correct, 
it not only provides a much-needed probe of the gas dynamics at 
intermediate radii in the galactic center, it also links the dynamics 
at event horizon-scales to larger radii approaching (and perhaps 
exceeding) the Bondi radius. 


APPENDIX A: INFERRING ANGULAR MOMENTUM 
AND TORQUE VECTORS FROM OBSERVATIONS 

The angular momentum and torque vectors used in section 2 depend 
simultaneously on /, Q, a, and ^ for G1 and G2, along with the 


time delay between their two orbits; quantifying these torques this 
requires knowledge of the correlated uncertainties in all of these 
parameters. We infer the torques directly from the data using an 
MCMC simulation. We use the emcee (Foreman-Mackey et al. 
2013) implementation of the affine-invariant algorithm presented in 
Goodman & Weare (2010). We tried both a standard statistic 
and a maximum-likelihood analysis in which the log-likelihood of 
a given orbit is proportional to 


InL = 


1 y Cvi - Xjf 

2 4 " +(/>’-)" 


+ \n[a-] + (fyif 


(Al) 


where the parameter / represents additional uncertainty in the data 
due to systematic errors; we find similar results in both cases. Our 
results are consistent with Gillessen et al. (2013b) and Pfuhl et al. 
(2015); table Al and figure A2 show a comparison. 

In order to convert from Kepler elements to Cartesian coor¬ 
dinates to compute likelihood functions, we define unit vectors 
along the eccentricity axis (pointing towards pericenter), along the 
angular momentum axis (orthogonal to the orbital plane), and the 
semi-minor axis (Murray & Dermott 1999). 


e - 


cos(du) cos(Q') - sin(du) cos(/') sin(Q') 
cos(du) sin(Q') -h sin(du) cos(/') cos(Q') 
sin(du) sin(/') 


7 = 


sin(/') sin(G') 
- sin(/') cos(Q') 
cos(/') 


b = 


- sin(du) cos(Q') - cos(du) cos(/') sin(Q') ' 

- sin(du) sin(Q') -h cos(du) cos(/') cos(Q') 

sin(/') cos(du) ^ 


(A2a) 

(A2b) 


(A2c) 


Note that this convention differs from that used by, e. g. Gillessen 
et al. (2012) and Phifer et al. (2013). We convert between the two 
coordinate systems with the following transformations: 


cos i' = jzl j = - cos i (A3a) 

(z X j)y 

tanfi' = ^^ = coin (A3b) 

We obtain positions and velocities using (Murray & Dermott 
1999): 


r = a(cos E - e)e + a Vl - sin Eb 


I - e cos E 


cos Eb - 


sin Ee ^, 


(A4a) 

(A4b) 


where the eccentric anomaly E is related to the mean anomaly via 
M - E - esinE, and the mean anomaly is related to time t via 
M = (t- fperi) dGM/a\ 

We do not fit for the uncertainty in the mass or velocity of 
Sgr A*, nor to the distance to the galactic center; instead, we adopt 
the best-fit values from Gillessen et al. (2009) and hold these con¬ 
stant throughout our analysis. 


ACKNOWLEDGMENTS 

We thank James Guillochon, Ramesh Narayan, and Eliot Quataert 
for interesting and helpful conversations which shaped the direc¬ 
tion of this work. Michael Johnson provided helpful comments, 
particularly relating to the EHT orientation measurement. A.-M.M. 
is supported by the National Aeronautics and Space Administra¬ 
tion through Einstein Postdoctoral Fellowship Award Number PF2- 
130095 issued by the Chandra x-ray Observatory Center, which is 

© 2015 RAS, MNRAS 000, 1-12 



Drag Force on Gas Clouds 11 


Table Al. Comparison of our derived Kepler elements with those published in Pfuhl et al. (2015) 
and Gillessen et al. (2013b). 


year (JD) 

«(”) 

e 

i' C) 

D' (°) 

a;(°) 

^peri 

G1 

0.5 + 0.1 

0.89 +0.03 

73 + 1 

22 + 5 

114 + 4 

2001.7 + 0.3 

gU 

0.4 + 0.2 

0.86 +0.05 

72 + 2 

21+5 

109 + 8 

2001.6 + 0.4 

G2 

1.2+ 0.2 

0.976 + 0.007 

62 + 2 

9 + 4 

94 + 2 

2014.22 + 0.08 

02^ 

1.0 ±0.2 

0.976 + 0.007 

62 + 2 

8 + 4 

97 + 2 

2014.25 + 0.06 


1 Pfuhl et al. 2014 
^ Gillessen et al. 2013b 
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Figure A1 - continued Keplerian fits to G1 and G2, with comparisons to 
astrometry and line-of-sight velocity measurements published in Gillessen 
et al. (2013b) and in Pfuhl et al. (2015). From top to bottom, rows show right 
ascension (RA), declination (dec), and line-of-sight velocity (vios) as func¬ 
tions of time. Below each plot, we show a residual after subtracting the best- 
fit model; this enables an easier comparison with the data. Blue regions show 
the 1- and 2-cr uncertainty in our fit to G2, and yellow curves show the same 
for Gl. The dark blue points show Br-y astrometry and velocity measure¬ 
ments for G2 published in Gillessen et al. (2013b). These are publicly avail¬ 
able here: https://wiki.mpe.mpg.de/gascloud/PlotsNData. The 
red points show measurements of Gl published in Pfuhl et al. (2015); 
these include both Br-y and L-band astrometry measurements. Light blue 
and orange points are the same as the dark blue and red ones, but show the 
larger error-bars determined in our maximum-likelihood analysis (via the / 
parameter in equation Al; see section 3.1). 
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Figure A2 - continued Probability distributions for Keplerian orbital ele¬ 
ments for G1 {top) and G2 (bottom). Each colored panel shows a projection 
of the seven-dimensional probability distribution onto the two-dimensional 
space defined by the indicated coordinates; at the end of each column, the 
green histogram shows the one-dimensional probability distribution of each 
variable. The correlated uncertainties between i' and Q' are essential for 
quantifying the uncertainty in j and thus t±_ in section 2. Similarly, the 
correlated uncertainties among all of a, e, i', and Q.' are needed to quantify 
the uncertainty in j and r. 
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